Band-center anomaly of the conductance distribution in one-dimensional Anderson localization 
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We analyze the conductance distribution function in the one-dimensional Anderson model of localization, for 
arbitrary energy. For energy at the band center the distribution function deviates from the universal form as- 
sumed in single-parameter scaling theory. A direct link to the break-down of the random-phase approximation is 
established. Our findings are confirmed by a parameter-free comparison to the results of numerical simulations. 
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The spatial localization of waves in a disordered potential 
can be considered as the most dramatic effect of multiple co- 
herent wave scattering ^ |[|. Due to systematic constructive 
interference in some part of the medium the wave function 
is spatially confined and decays exponentially as one moves 
away from the localization center 0|. The localization 
length Zioc can be probed non-invasively from the decay of the 
transmission coefficient (the dimensionless conductance [|]|) 
g, in terms of the average 



C\ = (-Ing) =2LAoc + 0(LO) 



(1) 



for system length L > Zioc Localization results in insulat- 
ing behavior of disordered solids at low temperatures [||, 4], 
and also can be realized in electromagnetic waveguides [7], 
where it is considered as an efficient feedback mechanism for 
lasing in disordered active media [^j]. 

One of the cornerstones of the theoretical understanding 
of localization is the universal approach of single-parameter 
scahng (SPS) [l^, [Tl|. In this theory it is assumed that the 
complete distribution function P{g) of the conductance can 
be parameterized by the single free parameter Ci . The de- 
pendence of Ci [and hence of P{g)] on L is then found from 
solving a scaling equation dCi/d{lnL) = P{Ci), where the 
universal scaling function /? does not depend on L, nor on 
any microscopic parameter (like the Fermi wavelength Ap, the 
transport mean free path Ztr, or the lattice constant a). 

The distribution function P{g) is completely determined by 
the cumulants 



a, ^ {{{- In gr)), 



(2) 



which are obtained as the expansion coefficients of the gener- 
ating function 



n=l 



n! 



(3) 



The first three cumulants are given by Eq. ([l|) for Ci, C2 = 
var In g, and C3 = (( (In g) — In g)^). The SPS hypothesis can 
then be phrased like this: All cumulants are universal func- 
tions ofCi. In the localized regime (Ci ^ 1), the universal 
SPS relations take the simple form | 



(4) 



These conditions are much more restrictive than the general 
upper bound C„ = 0(L/l\oc) from the theory of large- 
deviation statistics [ [l2| , [l3| ]: SPS assumes a lognormal distri- 
bution of g, with the variance of In g determined by the mean 
via the universal relation var In g = — 2(ln g). It is the viola- 
tion of this relation which frequently is used to indicated the 
break-down of SPS theory (see, e.g., Ref. [|l^ |l5||). 

In this paper we investigate P{g) in the most-studied 
and best-understood paradigm of localization, the one- 
dimensional Anderson model defined by the Schrodinger 
equation 



(5) 



on a linear chain of L sites (lattice constant a ~ 1) and a 
random potential with (V/) = and (ViVm) = 2DSi,n. The 
strength D of the potential fluctuations is taken to be small. 
We analytically calculate the cumulants Cn in the localized 
regime, with main focus on the energy region <C 1 around 
the band center of the disorder-free system. For E — we 
find the values 



C2/C1 = 2.094, C3/C1 = 0.568. 



(6) 



Cn/Ci^din+ 262n + OiL-^). 



The ratios C„/Ci with the higher cumulants also are finite. 
Hence P{g) complies with the restrictions of large-deviation 
statistics, but deviates from the special lognormal form as- 
sumed in SPS theory (this form is restored for \E\ > D). 

The conditions for validity of SPS have been a constant 
subject of intense debate [[uj, |l^ Originafly, SPS was 
derived within the random-phase approximation (RPA) for the 
scattering phase between consecutive scattering events [l^] . 
In the Anderson model the RPA is known to fail around the en- 
ergies E = ±2 (the band edges of the disorder-free system) 
|[r^], where Ap ^ ^tr- Indeed, the SPS relations (^ are vio- 
lated for all cumulants when one comes close to the band edge 
(2 — \E ^ D^/^) [p^, in coincidence with the expectations 

The RPA is also known to break down for the band-center 
case E — JT^]. However, the only consequence observed so 
far has been a weak anomaly in the energy-dependence of Zioc 
(hence, also of Ci) |^, ^], which differs at i? = by about 
9% from the predictions of perturbation theory [^. Surpris- 
ingly, the violation (^ of the SPS relations (Q) has not been 
noticed — quite the contrary, the relevance of the RPA for SPS 



2 



recently has been contested [p^ within an investigation of 
the Lloyd model, given by Eq. (g) with a Cauchy distribution 
for the potential | ]l6| , p3| ]. However, results obtained for the 
Lloyd model are not conclusive for the Anderson model and 
SPS, because in the Lloyd model formally D = oo and one 
encounters the modified universal relations C2 /Ci = 4 7^ 2, 
while Zioc varies smoothly with energy even around E = Q 
I p^ ] . Moreover, the higher cumulants have not been investi- 
gated. In previous numerical studies, the violations may have 
passed unnoticed because the small deviation of C2/C1 from 
the SPS value probably was not considered to be significant, 
and again the higher cumulants have not been investigated. In 
this paper, we also will establish a direct link between SPS 
and RPA. 

We now present the analytical calculation of the cumulants 
Cn of — In 5 in the vicinity of the band-center energy E = Q 
of the Anderson model, Eq. (|]). As pointed out many years 
ago by Borland [||], the dimensionless conductance g in the lo- 
calized regime is statistically equivalent to V'l'^' where tpL is 
the solution of the Schrodinger equation (|]) with generic ini- 
tial conditions tpo, tpi = 0{1). Because Ap ~ 4a, it is useful 
to introduce two slowly varying fields = 1)'/^ when 
I is even, x(0 = when / is odd, which can be 

considered as continuous functions with Langevin equations 



'/, = l(W + E^. (7) 



Here U and W independently fluctuate with (U) — 0, 
{U{Li)U{L2)) = ADSiLi - L2), and analogously for W. 

In order to calculate the wave-function decay and its fluctu- 
ations it is convenient to switch to the variables 



2x 



20 



(8) 



which are symmetric in and x- In the localized regime, u = 
— In (7 characterizes the global decay of the wave function, 
while the variable a (parameterizing the local fluctuations) is 
identical to the scattering phase of the reflection amplitude 
r — {ipL-i + i^^L)/ {ipL-i — iipL)- This parameterization 
allows us to draw a direct relation between SPS and RPA: SPS 
will turn out to be valid when a is uniformly distributed over 
(0,27r). 

The Langevin equations now can be translated into 
a Fokker-Planck equation for the joint distribution function 
P{u,a;x). For the sake of a compact presentation we use 
short-hand notations for the functions Sq, = sin a, Cq. = 
cos a, and introduce the rescaled position x = DL, as well 
as the rescaled energy e = E/D. The Fokker-Planck equa- 
tion then takes the form 

dxPiu,a]x) = 

i^l + du {sldu -cl + 2daSaCa) " edc.]P{u, a; x), (9) 
with the linear differential operator Ca — da{l + c^Y^"^ ■ 



The behavior of P{u, a; x) for large x can be analyzed by 
introducing into Eq. (ph the ansatz 



P{u, a; x) 



k=0 



T expinkiOx - £.u)fk{^, a), 

(10) 



where we require periodicity of fk {£,, a) in a. It then follows 
that the functions fk (C, ct) solve the eigenvalue equation 

flkfk^[Cl~eda.+acl-2da.Sa^Co.)+esl]fk, dD 

in which ^ appears as a parameter and fikiO is the fcth 
eigenvalue (arranged in descending order). In the vicinity of 
^ = 0, there is a finite gap between the largest eigenvalue no 
[which vanishes for ^ = 0, because of the normalization of 
P{u, a; x)] and fii. According to Eq. (|lO|), the asymptotic be- 
havior of the distribution function P{u, a; x) for large x hence 
is governed by fiQ, up to exponentially small corrections. A 
formal calculation of the moments of u (i.e., of — In 5) shows 
that the cumulant-generating function (^ is directly given by 
!]{£_) = x^oiO- Hence, 



C'n 



(12) 



where we expanded /io(C) — J2n=i M into a power se- 
ries. 

The expansion coefficients /i'"' can be calculated recur- 
sively for increasing order n by solving the hierarchy of equa- 
tions 



+ (13) 



fe=0 



which results when one introduces into Eq. (jllh the power 
expansions for ^0 and for fo{C,ct) = X^J^o /^^^^Q^)?"- I" 
each order n, we first integrate over a from to 2tt, which 
eliminates /^^"^ and hence gives /i^^"^ in terms of the quantities 
and /i^™) with m < n. Afterwards /^"^ can be obtained 
from Eq. ( [l3| ) by two integrations. The iteration is initiated 
for n = with /i'"^ = 0. This completely solves the problem 
to calculate the cumulants C„ in the localized regime. 

Let us illustrate the procedure for i? = 0. To start the itera- 
tion we consider Eq. ( pj| ) with n = 0, given by £^/*^°-' — 0. 
This differential equation is solved by the normalized function 



/27r 



r2(l/4)Vl + cos2a 



(14) 



which is identical to the stationary limiting-distribution func- 
tion lim2:^oo '^'^ P{(^^ x) of the variable a. 

Now the next iteration. Equation ( p3[ ) with n = 1 is given 

by 



(15) 
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FIG. 1: First three cumulants Cn = ((( — In g)")) for energy _E = 
in the Anderson model with D = 1/150, as function of system 
length L. The data points are the result of a numerical simulation. 
The slopes of the straight lines follow the predictions of Eq. (po|). 
The localization length li^c is taken from Eq. (p^). 



We first determine 



(1) 



27r 



r2(3/4) 



/o r2(i/4)- 

The prediction for the inverse locahzation length 

;ioc = r2(l/4)/[2OT2(3/4)], 



(16) 



(17) 



obtained by combining Eq. ( [l6| ) with Eqs. ^ and dl^), is 
identical to the result found in Refs. [ ^0| , |TJ. Then we solve 
for 



/W(a) = (l + c^)-V2 / dp{l + clr"' 



Jo 



From the next iteration n = 2 we obtain 



(18) 



da[(ci-A.W)/W(a) + 4/^°)(a)] (19) 







and also f^^'' (a). Analogously we 

obtain ^(3) withEq. ([l|), 
this is sufficient to determine the values for the first three cu- 
mulants 

Ci = 0.4569 DL, C2 = 0.9570 DL, C3 = 0.2595 DL. 

(20) 

They correspond to the anomalous ratios given in Eq. (||). 

The analysis of Eq. ( pj| ) can be straightforwardly caiTied 
out also for finite E/D. For E/D ^ 1, the stationary 
limiting-distribution function of a is given by /'•'^-'(a) — 
l/(27r), corresponding to a completely random phase. For 
n = 1 we find the coefficient /x'^' — 1/2, and the perturbative 
result Zioc = 4/Z3 is recovered 0]. In the next iteration we 
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FIG. 2: Same as Fig. [l], but for energy i5 = 0.1 (upper panel) and 
E — 2 (lower panel). The straight lines in the upper panel follow the 
predictions of perturbation theory [ p2| ] and single-parameter scaling 
I IC |. The straight lines in the lower panel are the predictions of Ref. 
1 17 1 (see text). 



obtain /i*^^^ = 1/2, while the higher coefficients all vanish. 
According to Eq. (|l2|), the SPS relations (^ then are reestab- 
lished. 

We have tested the predictions of the analytical theory 
against the result of a direct numerical computation of the con- 
ductance g for the Anderson model (||), by recursively increas- 
ing the length of the wire p4|]. The potential Vi was drawn 
independently for each site from a box distribution with uni- 
form probability 1 / V24D over the interval [— \/6D, V6D] . 
The data shown in the plots was obtained for D = 1/150 
(identical results are obtained for a Gaussian distribution with 
the same variance D). The cumulants were determined by 
averaging over 10^ disorder realizations. 

The result of this computation for the first three cumulants 
and £' = is shown in Fig. |l} The cumulants all increase 
linearly with the length L of the wire, and the slopes agree 
perfectly with Eq. ( pO[ ) [hence the localization length agrees 
with Eq. (|l7|) and the ratios of cumulants agree with Eq. (^]. 
The comparison is free of any adjustable parameter. 

For contrast, the upper panel of Fig. ^ shows the first three 
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tial. However, it can be questioned whether this universahty 
also extends to additional disorder in the hopping rates, since 
it is well known that the extreme case of purely off-diagonal 
disorder results in delocalization at E' = [ESl] . 



E 



FIG. 3: Energy dependence of the ratios of cumulants C2/C1 and 
C3/C1. The inset shows Ci in units of the perturbative resuU DL/2. 
The data points are the resuU of a numerical simulation of the An- 
derson model with D — 1/150. The curves are the analytical pre- 
dictions of this paper (E < 0.1), of perturbation theory [ p2| and 
single-parameter scaling (0.1 < E < 1.8), and of Ref. 
{E > 1.8). 



cumulants at energy E = 0.1, where the SPS relations 
hold and Ci = DL/2 follows from perturbation theory [E3]. 
The lower panel shows the results at the band edge E — 2, 
which are compared to the predictions Ci — 0.7295 D^/^L, 
C2 = 1.602 D^^^L, C3 = 0.7801 ZJi/^x of Ref. l^. 

In Fig. ^ we show the ratios of cumulants C2/C1 and 
C3/C1 as a function of energy. The inset shows Ci. The 
anomalous region extends up to i? ~ lOD. Around the band 
edge, the violations set in for 2 — i? < 3D^/^. Again, per- 
fect agreement is found between our analytical theory and the 
results of the numerical simulations. 

In summary, we have presented an analytical theory for the 
distribution function P of the dimensionless conductance g in 
the localized regime of the Anderson model, Eq. (^). The re- 
lations implied by single-parameter scaling theory for the 
cumulants C„ of — In g are violated not only around the band 
edges \E\ = 2, but also at the band-center energy E — 0, 
where the correct values are given by Eq. (^). Since the 
random-phase approximation is known to break down in both 
cases, our findings reestablish the relevance of this approxi- 
mation for single-parameter scaling, which recently has been 
contested QHll. 



Whether the single-parameter scaling hypothesis itself 
breaks down at E' = 0, or just persists in modified form, is 
an open question. The ratios (||) still imply universal relations 
between the cumulants for weak on-site disorder, i.e., they do 
not depend on the distribution function of the random poten- 
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